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On distribution modulo 1 of the sum of powers of a Salem 

number 

Dragan Stankov 


Abstract. Let 0 be a Salem number and P{x) a polynomial with integer 
coefficients. It is well-known that the sequence (9'^) modulo 1 is dense but 
not uniformly distributed. In this article we discuss the sequence (P(0”')) 
modulo 1. Our first approach is computational and consists in estimating the 
number of n so that the fractional part of {P(6'^)) falls into a subinterval of 
the partition of [0,1]. If Salem number is of degree 4 we can obtain explicit 
density function of the sequence, using an algorithm which is also given. Some 
examples confirm that these two approaches give the same result. 


1. Introduction 

Studying the distribution modulo 1 of the powers of a fixed real number 9 
greater than 1, has been of interest for some time. In his monograph [3, R. Salem 
considered the case of certain special real numbers 9. For instance, he showed that 
{6*"} tends to 0 in K/Z when 6* is a Pisot number. If 6* is a Salem number, {0"} 
is dense in M/Z, i.e. the fractional parts of 0” are dense in the interval [0, 1) but 
not uniformly distributed. (See [^, p. 87-89.) Moreover, Salem numbers are the 
only known numbers whose powers are dense in K/Z. Recall that a Pisot number 
is a real algebraic integer greater than 1 whose conjugates other than itself have 
modulus less than I. A Salem number is a real algebraic integer greater than 1 
whose conjugates other than itself have modulus less than or equal to I and at 
least one conjugate has modulus equal to 1. It is well known that one and only one 
of these conjugates 0“^ is inside the unit disc while the others are on the boundary. 
The degree 2t of 0 is necessarily even and at least equal to 4. 

We will use the following standard notation: 

Definition. For any real number x, we denote: 

(1) Integer part of x: [x] = max{n G Z : n < x\. 

(2) Fractional part of x: {x} = x — [x]. 

(3) Congruence modulo 1: x = x' mod 1 x — x' G Z. 

(4) Distance from x to the nearest integer: ||x|| = min{|x — n| : n G Z}. 

Let u = {u(n)) be an infinite sequence of real numbers. For 0 < x < 1, we 
define the repartition function 

Definition. 

/(x) = lim 2(;.a,rd{n < V|{M(n)} < xj 
N^oo N 


2010 Mathematics Subject Classification. Primary 11K06; Secondary 11R06. 


1 



2 


DRAGAN STANKOV 


There is an analogy between the repartition function and the distribution func¬ 
tion as well as between f'{x) and density function in the Probability theory. 

From now on, we suppose that 0 is a Salem number, u{n) = P{9"') where 
P{x) is a polynomial with integer coefficients. Denote conjugates of 6 by 9~^, 
exp(±2z7ra;i),..., exp(±2i7rwt_i). The sum of an algebraic integer and its conju¬ 
gates is an integer and therefore for all n € N, 

t-i 

d" -I- -I- 2 ^ cos 2Trnujj = 0 (mod 1) 

i=i 

so that the distribution of 0” (mod 1) is essentially that of —cos27rnWj. 

2. The main theorem 

The aim of this paper is to define an algorithm for determination of the repar¬ 
tition function and its first derivative, in the case 0 is a Salem number of degree 
four. The algorithm is described in the following 

Theorem 2.1. Let 9 be a Salem number of degree f, n, m positive inte¬ 
ger numbers, Pm(x) = ® polynomial of degree m with integer coeffi¬ 

cients. Let Qf^{x),Qf^{x),... ,Qff{x) be all branches of the inverse function of 
Qmix) = — 2 T^estricted to the domain [—1,1] where Tj{x) is Chebi- 

shev polynomial of the first kind of degree j. If domain of Q'f^{x) is [ak,l3k], 
Oik < Pk, k = 1,2,... ,K then we define its extension Sk(x) on M.' 

( Qk^ix); X e [ak,Pk] 

Qf:^{ak)-, X e {-oo,ak) , 

Qk^iPk); X e (/3fe,oo) 

in the case Q^^{x) is decreasing for x € [ak,Pk], 

{ -Qk^ix); X e [ak,Pk] 

-Qk^icik)] X e (-00,Q;fc) 

-Qk^iPk); X G i/3k,oo) 

in the case Q^"^{x) is increasing for x G [ak,l3k]- 

Let us denote g{x) = 7r“^ and let M be a natural number 

such that M > max({|Q;i|, |a 2 |j • ■ •, |<ta:|} Uil/^il) l/^zj) • • •, |/3 k'|})- The repartition 
function of u{n) = is f{x) = Yh^-m {g{x -I- *) — g{i)). The first derivative 

of the repartition function is 

M K 

(2.1) f{x)=n-^ ^|(arccos(Qfc\T-f t)))'|. 

i— — M k—1 

Proof. Let conjugates of 9 be 9~^, exp(2t7rw), exp(—2i7ru;). Since for any 
natural n 

Oj (0"-^ -f 0“"^ -I- 2 cos 2TTnjuj) = 0 (mod 1) j = 0,1,..., m 

m m 

Pm{d^) = 'Yj ~ ^ 2Trnjuj) (mod 1) 

j=o j=o 

so that the distribution of Pm{9'^) (mod 1) is essentially that of 

m m 

—2 Y^ o,j cos 2'Knju} = —2 OjTj (cos 2Tmuj) = Qm (cos 2'Knuj), 

1=0 j=o 



ON DISTRIBUTION MODULO 1 OF THE SUM OF POWERS OF A SALEM NUMBER 3 


where Tj{x) = J2k=o^k'^ Chebishev polynomial of the first kind. Hence we 
have 


Qm{w) = = -2(cmu;™ + Cm-iw" 

j—0 k—0 

where we denoted w = cos 2Trnuj and 

„ _ „ him) 

— ^m^rn ’ 

,(m—1) . A'm) 

Cm—1 — ^m—l^rn—1 


Co). 


( 2 . 2 ) 


Co = a 




+ • • ■ + O-mb 


(m) 


Co — ^ 0 ^ 0 ^^ + 0,1^0 ^ + 02^0 ^ + . • . + am^o ^ 

If we denote an integer M = 2 jyjLo |cj I then it is obvious that Qm{w) G [—M, M] 
thus {Qm{w)} < X there is an i G {—M, —M+1,... ,M} such that i < Qm{w) < 
i + X. Now we conclude that there is a /c G {1, 2,..., K} such that Q^^{i) < w < 
Q^^{i + x) or > w > Q^^(i + x). Previous two inequalities can be replaced 

with Sk{i) > w > Sk{i + x). It is easy to verify that (27r)“^ arccos(cos27rt) = ||t|| 
and, as a consequence of this, that (27r)“^ arccos(— cos 27rt) = (27r)“^ arccos cos(27rt+ 
tt) = \\t + 1/2||, for that reason we have 

(27r)“^ arccos((5fc ^(i)) < ||na;|| < (27r)“^ arccos(Q^^(i + t)) 
in the case ^(*) is decreasing, or 

(27r)“^ arccos(—Q^^(i)) < ||na; + I/2|| < (27r)“^ arccos(—+ x)) 
in the case is increasing. 

It is fulfilled that ||na;|| and ||nw + l/2|| are uniformly distributed on [0,1/2] be¬ 
cause I, uj are Q-linearly independent [2], Theorem 5.3.2 so we can use [2] Theorem 
4.6.3. Consequently, for all L, R such that 0 < L < i? < 1/2, 

lim ^cardjn < NIL < ||na;|| < i?} = 2{R — L), 

N—^oo N 

lim ^cardjn < N\L < \\noj -f ^|| < R} = 2{R — L). 

N —>00 N 2 

Let Ki < K he natural number such that <5/^, <5^^,. •., are decreasing and 
Q~Ki+i^ Qifi+ 2 J • ■ • j Qli increasing. Now we can determine the repartition func¬ 
tion ^ 

fix) = lim —cardIn < iV|{(5m(cos(27rna;))} < t} = 

N^oo N 

1 M K^ 

TV™ U (U{- < iV|(27r) ^ arccos((5^ ^(i)) < 

i— — M k—1 

< \\nuj\\ < (27r)“^ arccos(Q^^(z + x))} [J 

K 

[J {n < A^|(27r)“^ arccos(—< 
k^Kx + l 

< \\nu)-\-l/2\\ < (27r)“^ arccos(-Q^^(z + x))} ^ = 

M K 

E E (tt ^ arccos(S'fc(i-f a;)) — tt ^ arccos)^^)^))) = 

i— — M k—1 
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M 

(gix + i) - 9{i)) 

i^-M 

because all sets in the double union are disjoint. Now it is obvious that the first 
derivative of the repartition function is 

M 

(2.3) f'{x)= Y 9 '{x + i). 

i=-M 

□ 

Remark. Since {P{d^)} = {P(0") + Z}, Z £ Z we can take, without loss of 
generality, that oq = 0. 

Remark. Since g{x + i) — 9 ( 1 ) = 9 {x + i) + c — g{i) — c it is clear that f{x) 
will not be changed if we take g{x) + c, c G K instead of g{x). 

Remark. It is necessary to introduce branches of the inverse function Q{x) 
more precisely. Since Q{x) is a polynomial we can introduce a partition of [—1,1] 
— 1 = a;o < < • • • < xk = 1 such that Q'{xi) = Q'{x 2 ) = ■ ■ ■ = Q'{xk-i) = 0 

and Q'{x) is positive or negative on each sub-interval {xk-i,Xk), k = 1,2,..., K. 
Let us introduce ak,f3k' if Q'{x) is positive on (xk-i,Xk) then Q{xk-i) = a*, 
Q{xk+i) = Pk', if Q'ix) is negative on {xk-i,Xk) then Q{xk-i) = Pk, Q{xk) = a/c- 
Now we define Qk^{x) as the inverse function of Q{x) on [xk-i,Xk\- Let us notice 
that the first derivative of Qk^{x) can tends to infinity only in end points of its 
domain [ak^Pk]- 

Corollary 2.1. Let the line x = v, v € [0,1] be a vertical asymptote of the 
graph of the first derivative of the repartition function y = f{x). lima;_>^_o fix) = 
00 if and only if v = 1 or v = {Pk}- lima;_,.^+o fix) = 00 if and only if v = 0 or 
V = {ofe}, k=l,2,...,K. 

Proof. We proved in the Theorem |2.1| that 

M M K 

fix)= 'Y 9 'i^+ ’>■)= 'Y aiccos'jSkji + x)) = 

i——M i— — M k—1 

M K 

= -7r-^ Y Y^^~ ^kii + x))-^/'^S'kii + x) = 

i— — M k—1 
M K 

= -7r-i Y + 

i— — M k—1 

lim 2 ;_>„_o fix) = 00 if and only if there are iq, ka such that 
- iQla)'^if+x))~^^'^S'kfio + x) = 00 . 

There are two cases: either Qk^if + f) = ±1 or lim 2 ,_).„_o <5'^^ (*o + x) = 00 . If 
Q^^iio + v) = ±1 then iQ — v = Qkgi^l)- Since (^^^(±1) is an integer and v G [0,1] 
we conclude that either n = 0 or w = 1. But = 0 is impossible because z; — 0 
will be out of the domain of Qfi^ix). If lim 2 ,_,.i,_o S^^^iio -|- cc) = 00 then, using the 
last remark, either zq -I- = akg or ig + v = Pko ■ Again + v = akg is impossible 

because Zq + z; — 0 will be out of the domain of S'f,^ix). If zg -I- z; = Pko then we 
conclude that v = {Pko} with an exception: if Pko is an integer then its fractional 
part is 0 but, as we have seen, z; = 0 is impossible. Nevertheless the claim is true 
because, in that case the line x = 1 should be a vertical asymptote of the graph. 

It is obvious that lima;_,.„_|_o fix) = 00 if and only if z; = 0 or z; = {as,} can be 
proved completely analogously. □ 
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We present next procedure for sketching the graph of the first derivative of the 
repartition function, resulting from the previous corollary: 

(1) Find set A of local minimum points, set B of local maximum points and 
set S of (horizontal inflection) stationary points of Q{cos{t)). 

(2) Find set A of fractional parts of values at local minimum points, set B of 
fractional parts of values at local maximum points and set S of fractional 
parts of values at stationary points of Q(cos(t)). 

(3) Let xo, Xi,..., Xr, 0 = xq < Xi < ■ ■ ■ < Xr = i he sorted r + 1 elements of 
AUBUS. 

(4) If e >1 U 5*, Xij^i € BVJ S then f'{x) has vertical asymptotes x = Xi, 
X = Xi+i on interval {xi,Xi+i) so f (x) has the shape of U. 

(5) If Xi G AUS”, Xi+x G A\{BLiS) then f'{x) has vertical asymptote x = Xi 
on interval {xi, x^+i) so f’{x) has the shape of left half of U, we will denote 
it by [. 

(6) If Xi G B \ {A LI S), cci+i G B L S then f'{x) has vertical asymptote 
X = Xi+i on interval {xi,Xi+i) so f{x) has the shape of right half of U, 
we will denote it by J. 

(7) If iCi G B\(Ails'), Xi+i G A\(BUS) then/'( cc) has no vertical asymptote 
on interval (xi^Xi^-i) so f'{x) has the shape of 

Remark. In the first item of the previous procedure we have to solve the 
equation —Q'(cost) sint = 0 sint = 0 V Q'(cost) = 0. Solutions of sint = 0 
are t = kir, k G Z. If a; is a solution of Q'{x) = 0 and —I < x < 1 then t = 
arccos(x) + 2fc7r, /c G Z is a stationary point of Q{cos{t)). Thus, if a solution of 
Q'{x) = 0 is out of K or greater than I in modulus we should ignore it. 

In the second item we have to find {(^(cosfcTr)} = {(5(±1)} = 0, fc G Z so that 
0 G A U B U S. Since the fractional part of a real number is in [0,1) we should take 
that 0 and I must be both in or both out of set A, as well as B and S. We conclude 
that I G A U B U B. 


3. Linear, quadratic and cubic polynomial 


If P{x) = aix then, using the notation of the Theorem 2.1 we have only one 
branch of the inverse function of Q{x) = —2aix i.e. Q~^{x) = ^, so that g(x) = 

7r“^ arccos(^). The repartition function is f(x) = (arccos(—|^) — 

arccos(— 2 ^))- The first derivative of the repartition function is 


f'(^) = 


2ai7r 


I - 


(x—2ai 

4a? 


(x—2ai + l)^ 
^ 4a? 


+ ••• + 


(a;+2ai-l)2 

^ 4a? 


If we take oi = 1 we get Dupain’s formulae cited in 

Hereafte r we suppose that P(x) = a 2 X^ + aix and then, using the notation of 
we have two branches of the inverse function of Q(x) = —4a2X^ — 


the Theorem 


2.1 


2aix + 2«2 i.e. Q^\x) = ^ 

denote A — af + Sa^ — 4a2x; 

„ , I _ / ~ arccos((5)"^(x)), if A > 0 and — 1 < Qi^(x) < I, 

^ * ( 0, othervise; 


G2(x) = 


— arccos((52 ^(x)), if A > 0 and — I < Q 2 ^(^) — 1 

0, othervise; 


then we can prove that g(x) = tt ^(Gi(x) + G 2 (x) — tt) in the case that x < 
202 + 2ai — 402 and 02 > 0; g(x) = 7r“^(Gi(x) + G 2 (x) + tt) in the case that 
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X > 202 + 2ai — 402 and 02 < 0. Since 

f __ if A>0, |Q^^(a;)|<l 

= < (16a| —(ai+(a^+8a| —4o2a:)i/2)2)(a^+8a| —4a2a^)^/^ ’ 

1 0 , other vise; 

f _^ if A>0, |Q;'h 2 :)l<l 

= s (16a| —(ai —( 0 ^+ 803 —4a2a:)i/2)2)(a^+8a| —4a2a:)i/2 ’ 

1 0 , other vise; 

we can determine g'(a;) explicitly by g'{x) = ■k~^[G’i{x) + G’ 2 {x)). The repartition 
function is /(x) = + *) “ 5(*))) ^ = 4 |o 2 | + 2|oi| + 2 |a 2 | . The first 

derivative of the repartition function is f'{x) = d '+ 0 - 

Since = 0 Q{x) has extremum V = + 2 o 2 at x = —Using 

previous Corollary we can conclude that the graph of the first derivative of the 
repartition function y = f'{x) has an inner vertical asymptote x = v, v = {V} S 
( 0 , 1 ) if and only if 

(3.1) oi^O, U^Z. 

4a2 

If 02 > 0 then Q{x) has maximum V so that lim 2 ;_>.„_o/^(x) = 00 . In that case 



branches of the inverse function of Q{x) = —Soax^ — 4 a 2 X^ + ( 6 a 3 —2oi)x + 2 o 2 . Its 
explicit formulas are clumsy so we will not cite them here. The roots of Q'{x) = 0 
are 

—02 ± \/9a§ + 02 — SoiOa 

603 

In the Table using the procedure for sketching the graph of /'(x), we represent 
different shapes of graphs. 

Table 1. Procedure for sketching the graph of the first derivative 
of the repartition function 


03 , 02,01 

Xi 

X2 

Q{xi) 

Q{^i) 

A 

B 

S 

nx) 

1 , 1,1 

-0.61 

0.27 

- 0.11 

2.63 

.89,0,1 

.63,0,1 


U..U 

3,5,6 

- 0.68 

0.12 

4.22 

10.39 

. 22 , 0,1 

.39,0,1 


LUJ 

3,3,10 

-0.17 

-0.17 

6.11 

6.11 

0,1 

0,1 

.11 

uu 

1 ,- 1,-2 

-0.5 

0.83 

-5 

4.48 

0,1 

.48,0,1 


UJ 

1,2,3 

-0.67 

0 

2.82 

4 

.82,0,1 

0,1 


LU 

1 ,- 2,-2 

-0.39 

1.06 

- 6.21 


.79 

0,1 



1 , 2,-2 

-1.06 

0.39 


6.21 

0,1 

.21 



1 , 0,0 

-0.5 

0.5 

-2 

2 

0,1 

0,1 


U 

1,1,4 

^ K 

i M 



0,1 

0,1 


U 
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4. Some polynomials of degree m > 3 

If P{x) = a:™ then u{n) = (0™)”. We need a Salem’s lemma which is proved 

in mm- 

Lemma 4.1. If 6 is a Salem number of degree 2t then so is 9^ for all m gN. 

Thus, if we denote 9i = 0™ we get u{n) = 0” which is a well known sequence. 
Doche, Mendes France and Ruch [4] proved next 


Lemma 4.2. Let 6i be a Salem number of degree 2t, then the repartition function 
f{x) of the sequence ( 0 ") modulo 1 satisfies 


/'(a;) = 1 + 2 ^ Jo(4fc7r)‘ ^cos27rfcx 


k=l 


on ( 0 , 1 ) for all t > 2. 


Here Jo{-) is the Bessel function of the first kind of index 0. 

There are integer coefficients aj of P{x) such that Q{x) = —2"’'a;™. We can 
use equations 2.2 to find such aj. The power a;" can be expressed in terms of the 
Chebyshev polynomials of degrees up to n (for proof see [6], chapter 2.3.1): 

[m/2] 

(4.1) a:"* = 2^-’" 


k—0 




where the dash (^') denotes that the fc-th term in the sum is to be halved if m is 
even and k = m/2. Let m be odd, we conclude that if 


(m-l)/2 

p{x) = 


'i—2k 


then Q{x) = —2"*a:'", its inverse function can be easily found: Q ^{x) = — ’f/x/2. 
Using [ 2 T] we obtain that 


2™-l 


Va^ + * 


f(.x) = - Y — 

^ i=- 2 ™ 2m{x + — ( "^x + z)2/4 

We will show that /'(1/2 + x) = /'(1/2 — x), |x| < 1/2. For that reason x = 1/2 is 
a line of symmetry of the graph of f'{x). It is convenient to introduce 


g'ix) = 


1 


TT 2mxy/l — ( y/x)^/A’ 

then we have 

/'( 1/2 + x) = + a: + i) 

= Er-”2™ff'(-l/2-a:-*) 


(g'(x) is even) 
(* = J - 1 ) 


= Ei=-2™+i - X - j + 1) (j = -i) 

= Er-”2™ff'(l/2-a: + *) 

= /'(l/2-x). 


Let m be even, we conclude from |4.1| that if 

ml2-l 


p{x) = Y 




L — 2k 


+ 


1 / m 


2 \m/2 
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then Q{x) = —its inverse function can be easily found: Q = ± '^12. 
Using the algorithm presented in the Theorem |2.1| we obtain that 

2 ’" 


'TT 


V—a: + i 


i=i m(—x + i)^l — ( V—a; + iy/A 



0 OJ 0.4 0.6 0.8 


— .x,0.63I.0.887.x 
150 ■ 

Figure 1. f'{x) (the red curve) of ({^(0")}), where P{x) = + 

x"^ + X and distribution histogram of the sequence. The shape of 

fix) isU..U- 



— .ii,0Jl8.0.391,i 
152' 


Figure 2. f'{x) (the red curve) of ({P(d")}), where P{x) = 3a:^ + 
5a:^ + Qx and distribution histogram of the sequence. The shape of 

fix) is LUJ- 


5. Some examples 

To illustrate the algorithm and the procedure, we give examples of distributions 
for some polynomials from the Table 1. Using the definition of f{x) and the well 
known fact that the first derivative of the function f{x) on a small interval [a, 6] 
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Figure 3. f'{x) (the red curve) of ({P(0"')}), where P{x) = x^ — 
x^ — 2x and 0 is a Salem number of degree 4. The distribution 
histogram of i{P{df)}) where 0i is a Salem number of degree 6. 

We can notice that it is close to the density function of the uniform 
distribution. The shape of f'{x) is UJ. 

could be estimated with the finite difference: f{b) — f{a) divided by & — a, we will 
approximate f'{x). The interval [0,1] is divided into p pieces. We compute the 
fractional part of P(6>") for 1 < n < N, and count the number of n so that the 
fractional part of P(0”) falls into each of subintervals. The vertical axis indicates 
the number of such n divided by N/p so that this normalisation results in a relative 
histogram that is most similar to f'{x). 

Y. Dupain and J. Lesca [5] conclude that for large degrees t, the sequence 
({0”}) is close to being equidistributed, a fact that S. Akiyama and Y. Tanigawa 
[T] make very explicit in their article. We can notice that the same property is valid 
for the sequence {P{6^)) modulo 1 (Fig. 3). 
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